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The coexistence of ferroelectricity and magnetism in ABO3 perovskite oxides is rare, a phe¬ 
nomenon that has become known as the ferroelectric u d° rule”. Recently, the perovskite BiCoC>3 
has been shown experimentally to be isostructural with PbTiOs, while simultaneously the d 6 Co 3+ 
ion has a high spin ground state with C-type antiferromagnetic ordering. It has been suggested that 
the hybridization of Bi 6s states with the O 2 p valence band stabilizes the polar phase, however 
we have recently demonstrated that Co 3+ ions in the perovskite structure can facilitate a ferroelec¬ 
tric distortion via the Co 3d - O 2p covalent interaction (Weston et al. , Phys. Rev. Lett. 114, 
247601 (2015)). In this report, using accurate hybrid density functional calculations, we investigate 
the atomic, electronic and magnetic structure of BiCoOs to elucidate the origin of the multifer¬ 
roic state. To begin with, we perform a more general first-principles investigation of the role of d 
electrons in affecting the tendency for perovskite materials to exhibit a ferroelectric distortion; this 
is achieved via a qualitative trend study in artificial cubic and tetragonal LaBOs perovskites. We 
choose La as the A-cation so as to remove the effects of Bi 6s hybridization. The lattice instability 
is identified by the softening of phonon modes in the cubic phase, as well as by the energy lowering 
associated with a ferroelectric distortion. For the LaBOs series, where B is a d° — d 8 cation from 
the 3d-block, the trend study reveals that increasing the d orbital occupation initially removes the 
tendency for a polar distortion, as expected. However, for high spin d 5 — d 7 and d 8 cations a strong 
ferroelectric instability is recovered. This effect is explained in terms of increased pseudo Jahn-Teller 
(PJT) p — d vibronic coupling. The PJT effect is described by the competition between a stabilizing 
force (K 0 ) that favours the cubic phase, and a vibronic term that drives the ferroelectric state ( K v ). 
The recovery of the lattice instability for high spin d 5 — d 7 and d 8 cations is due to (i) a reduction 
in Ko due to a significant volume increase arising from population of the cr-bonded axial d e g or¬ 
bitals, and (ii) an increase in the K v contribution arising from increased p — d hybridization; our 
calculations suggest that the former mechanism is dominant. Surprisingly, we are able to show that, 
in some cases unpaired electron spins actually drive ferroelectricity , rather than inhibit it, which 
represents a shift in the understanding of how ferroelectricity and magnetism interact in perovskite 
oxides. It follows, that for the case of BiCoOs, the Co 3+ ion plays a major role in the ferroelectric 
lattice instability. Importantly, the ferroelectric polarization is greatly enhanced when the Co 3+ 
ion is in the high spin state, when compared to the nonmagnetic, low spin state, and a large cou¬ 
pling of the electric and magnetic polarization is present. Generally, for d 5 — d 7 B-cations in ABO3 
perovskites, an inherent and remarkably strong magnetoelectric coupling exists via the multiferroic 
crossover effect, whereby switching the spin state strongly affects the ferroelectric polarization, and 
potentially, manipulation of the polarization with an externally applied electric held could induce 
a spin state transition. This novel effect is demonstrated for BiCoOs, for which the ground spin 
state is switched by reducing the internal ferroelectric polarization. These results provide a deeper 
insight into perovskite ferroelectrics and multiferroics. 

PACS numbers: 75.50.Pp, 77.80.-e, 75.85.+t 


I. INTRODUCTION 

The apparent mutual exclusion of ferroelectricity and 
magnetism in ABO3 perovskite oxides has both fasci¬ 
nated and puzzled researchers in the held of multiferroics 
for decades BB- It is generally known that the presence 
of ^-electrons at the B-site cation, which is a requirement 
for magnetism, reduces or removes the tendency for non- 
centrosymmetric distortions of the B06 octahedra, a phe¬ 
nomenon which has been dubbed the u d° rule”. While 
there has been much progress in this held, the observed 
mutual exclusion between ferroelectricity and magnetism 
remains to be fully understood [ 3 !. 


Perovskite oxides undergo various symmetry lowering 
distortions which decrease the total energy of the system. 
In the case of a ferroelectric material, the noncentrosym- 
metric distortions of B06 octahedra leads to breaking of 
inversion symmetry and a net macroscopic electric po¬ 
larization [ 4 ]. It is more common, however, for ABO3 
perovskites to undergo centrosymmetric distortions in¬ 
volving the tilting and rotation of BC>6 octahedra [ 5 ]. 
It has been suggested that competition with such cen¬ 
trosymmetric distortions precipitates the ferroelectric d° 
rule [I]. 

Ferroelectricity in perovskite oxides can have differ¬ 
ent origins, e.g., from canted spins in frustrated magnets 
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[ 6 ], or from size and electrostatic effects in the hexago¬ 
nal phase [7]; however, in this work, by using the term 
ferroelectricity we are implicitly referring to “proper” fer- 
roelectricity. Here, p — d covalent bonding between the 
B-cation d states and surrounding O 2 p states can sta¬ 
bilize the ferroelectric phase, which is characterized by a 
shift of the B-cation sublattice with respect to that of the 
surrounding O ions |4j. Examples of proper ferroelectrics 
include the classic cases of BaTiOs and PbTiOs, where 
the Ti 4+ ion does indeed have a d° electron configura¬ 
tion, and therefore both BaTiOs and PbTiOs are non¬ 
magnetic ferroelectric perovskites, consistent with the d° 
rule [Tj. It has recently, however, been pointed out by 
Bersuker [ 8 ], that certain 3 d n (n > 0) magnetic cations 
should be unstable against a pseudo Jahn-Teller (PJT) 
driven ferroelectric distortion. This proposition suggests 
that there exists no fundamental explanation for the fer¬ 
roelectric d° rule within the PJT theory for ferroelec¬ 
tricity, although an examination of experimental results 
suggests, at least empirically, that magnetic B-cations 
(almost) never exhibit proper ferroelectricity. Neverthe¬ 
less, currently there is a strong research effort focused on 
the design of magnetoelectric multiferroics, as the cou¬ 
pling of ferroelectricity and magnetism within a single 
phase material promises novel and exciting opportunities 
for future spintronic devices mm- 

The perovskite oxide BiCoOs has been shown to ex¬ 
hibit a tetragonal ferroelectric phase as the ground state 
(i.e., isostructural with PbTiOs) fT5| . At the same time, 
the d 6 Co 3 + ion is in the high spin state, with C-type 
antiferromagnetic ordering. This system is therefore a 
rare example of a d n perovskite exhibiting proper ferro¬ 
electricity, and represents an opportunity to to investi¬ 
gate the interaction of ferroelectricity and magnetism in 
perovskite oxides. It has been proposed, based on first- 
principles calculations [Id, that the ferroelectric phase 
of BiCoOs is driven by hybridization of the Bi 6s states 
and the O 2 p valence band; however, we recently demon¬ 
strated that Co 3+ ions doped into PbTiOs exhibit novel 
ferroelectric and multiferroic effects H3, and we therefore 
propose that the Co 3+ ion is responsible for the strong 
ferroelectric distortions exhibited by BiCoOs HE as well 
as for the multiferroic effects predicted theoretically [181 • 


At this point, it is clear that more work is needed, both 
experimental and theoretical, so as to understand novel 
multiferroic materials such as BiCoOs, as well as to fur¬ 
ther elucidate the complex interactions of ferroelectricity 
and magnetism on a microscopic scale. In this report, the 
interaction of ferroelectricity and magnetism in ABO 3 
perovskites is investigated in detail. It is demonstrated 
that the presence of d-electrons at the B-cation does not 
necessarily rule out a strong ferroelectric distortion; in 
fact, in the case of BiCoOs, the presence of unpaired d- 
electrons actually drives ferroelectricity, rather than in¬ 
hibits it. The manuscript is presented as follows: the 
methodology is described in section|H|and our results and 


discussion are presented in section |III| In section III A 


we provide a brief introduction into perovskite ferroelec¬ 


tricity and multiferroicity. Section IIIB investigates the 
interaction of ferroelectricity and magnetism via a qual¬ 
itative trend study for artificial LaBOs perovskites, and 
in section ra we describe quantitatively the origins of 
ferroelectricity and multiferroicity for the realistic per¬ 
ovskite BiCoOs. Our key findings are summarized in 
section m 


II. METHODOLOGY 

Calculations are performed within the density func¬ 
tional theory (DFT), using the screened hybrid func¬ 
tional of Heyd, Scuseria and Ernzerhof (HSE) [19] for the 
exchange-correlation potential. We use predominantly 
the standard implementation HSE06, where the screen¬ 
ing length and mixing parameter are fixed to 10 A and 
0.25 respectively. While the inclusion of 25% Hartree- 
Fock in the calculation of the exchange potential gives 
a highly accurate account of the atomic and electronic 
structure [20;. for the calculations in Fig. [lOj the mix¬ 
ing parameter is reduced to 0.125 so as to give a more 
accurate description of the spin splitting energy mm- 
The valence electrons are separated from the core by use 
of the projector augmented wave (PAW) [23] pseudopo¬ 
tentials as implemented in the VASP package [24]. The 
energy cutoff for the plane wave basis set is 500 eV. For 
the 5-atom LaBOs primitive cells a 6x6x6 k -point grid 
is used for integrations over the Brillouin zone. For the 
20 -atom BiCoOs cells a 4x4x3 grid is used. Full opti¬ 
mization of the lattice vectors is allowed, and the internal 
coordinates are relaxed until the forces are less than 0.001 
eV/A. The spin state of each perovskite is controlled us¬ 
ing fixed-spin moment calculations. The phonon frequen¬ 
cies are calculated using the method of finite differences. 
As has been previously reported for phonon calculations 
with DFT-VASP [25] . we find that large values for the 
displacements are required to obtain reasonable results, 
and we use displacements of 0.04 in lattice units; dis¬ 
placements this large are expected to be outside of the 
linear regime, however we also present the results of total 
energy calculations and lattice relaxations to confirm the 
findings of the phonon study. 


III. RESULTS AND DISCUSSION 

A. Perovskite Ferroelectricity and the Pseudo 
Jahn-Teller Effect 

To begin with, we describe the origins of perovskite 
ferroelectricity in terms of the PJT theory. The cubic 
crystal structure of an ABO 3 perovskite oxide is shown 
in Fig.JlJa). The A-cation is 12-fold coordinated to oxy¬ 
gen, and generally, the interaction with the surrounding 
O ions is ionic, in that the A-cation valence electrons 
are transferred to the O ion 2 p band without significant 
hybridization. On the other hand, the B-cation (shown 
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FIG. 1: (Color online) (a) The unit cell of an ABO 3 per- 
ovskite oxide with a cubic crystal structure. The B-cation 
is at the centre of the unit cell and is 6 -fold coordinated to 
oxygen, (b) The same system with a tetragonal ferroelectric 
distortion. The distortion is characterized by an elongation of 
the c-axis compared to the in-plane lattice vectors, as well as 
a shift of the B-cation with respect to the surrounding O-ion 
octahedron; the internal distortions break inversion symme¬ 
try, leading to a net macroscopic electrical polarization. The 
A-cation of the perovskite structure is represented by a green 
sphere, the B-cation is blue, and O ions are red. 


here at the center of the unit cell) is 6-fold coordinated 
to oxygen, and the B-cation d states will generally show 
a significant covalent interaction with the O 2 p states. 
A basic molecular orbital diagram for this covalent in¬ 
teraction is shown in Fig. [2] The O 2 p and B-cation 
d states will form bonding and antibonding states, with 
the former having predominant O 2 p character, and the 
latter having predominant B-cation d character. Due to 
the octahedral crystal field, the free ion d band is split 
into a t 2 g triplet and an e g doublet, for which the p — d 
interaction is 7r-bonding and cr-bonding, respectively. 

Perovskite oxides can deviate from the ideal cubic ge¬ 
ometry presented in Fig. [lja). For the case of a proper 
ferroelectric distortion, the B-cation at the center of the 
cubic unit cell will shift with respect to the surrounding 
O-ion octahedron, leading to a net electrical polarization, 
as shown in Fig.[ljb). The driving force for the lattice in¬ 
stability is the PJT effect p26ti28l . and this phenomenon 
can lead to coupled magnetic and dielectric polarizations 
[HI ESI EH; here, we present the key points of the theory. 
For a perovskite oxide originally in the cubic (paraelec- 
tric) phase, for a symmetry-lowering ferroelectric distor¬ 
tion, the adiabatic potential energy surface E(q) along 
the ionic displacement q has a curvature K, 


K = 



(i) 


The lattice exhibits an instability if the curvature at 
the high symmetry point is negative, K < 0. As 
E = (T 0 | H | T 0 ), where H is the Hamiltonian and T 0 
is the corresponding wavefunction in the ground state, 
Eq-0 can be expanded, 



where, Tq = • From Eq. j^J it is clear that K is 

the sum of two competing terms, 


K = K 0 + K v 


where, 


( 3 ) 


K 0 = 




( 4 ) 


Within second order pertubation theory, K v is expressed 
as follows EH ESI, 


k„ = 


E 


IH(^ 


s «Jt 


T 

^ r 


(Eq — E n ) 


( 5 ) 


Here, T 0 is the ground state wavefunction of the unper¬ 
turbed Hamiltonian with energy Eq, and T n are excited 
states with energy E n . The PJT effect provides the driv¬ 
ing force for ferroelectric distortions in perovskite oxides 
[8]. Whether or not a polar ferroelectric phase is favoured 
(i.e., a negative curvature K ) depends on the relative 
strengths of the contributions Kq and K v in Eq. [3] The 
term Kq described by Eq. [4] is always positive - as the 
wavefunction T 0 represents the ground electronic state 
for the unperturbed Hamiltonian, moving the coordi¬ 
nates along q will always increase the total energy. There¬ 
fore Kq provides a restoring force that favours the undis¬ 
torted phase. The term K v is the so-called “vibronic 
term”, which lowers the total energy as it corresponds to 
the response of the electronic wavefunction to the distor¬ 
tion, which is the mixing of To with appropriate excited 
states, T n . The matrix elements (T 0 |(JiF/Jg)o|T n ) in 
Eq. [5] are only non-zero if To and T n have the same 
spin, and due to the denominator (Eq — E n ) it is clear 
that a strong polar distortion requires that low-lying ex¬ 
cited states exist, with the same spin-multiplicity as the 
ground state. For PJT driven ferroelectricity, one can 
envision three possible scenarios: (i) when Kq K v , 
the system will be stable against a ferroelectric distor¬ 
tion, (ii) when Kq ~ K v , the two phases have a similar 
energy, and (iii) when K v <C K v , the system will exhibit 
a strong lattice instability that favours the polar phase. 

For the case of perovskite oxides, the vibronic term of 
Eq. [5] describes how a ferroelectric distortion is driven by 
increased overlap of O 2 p states with the d states of the B- 
cation mm- From Eq. |5j it is clear that the strength of 
the vibronic coupling between a p and d orbital depends 
on, firstly, how strongly the overlap of the orbitals is 
increased by the distortion, (p\(SH/Sq)\d), and secondly, 
the energy gap between the two states (E p — Ed). In 
principle, the O 2 p states can form covalent bonds with 
the non-axial t^g orbitals, or the axial e g orbitals. The 
t 2 g states lie at a lower energy than the e g states, and 
so according to the denominator in Eq. [51 should lead 
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B-cation d 
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Perovskite 

Crystal structure 

Reference 

LaScOs 

Orthorhombic 

m 

LaTiOs 

Orthorhombic 

m 

LaV0 3 

Monoclinic 

m 

LaCrOs 

Orthorhombic 

m 

LaMnOs 

Orthorhombic 

EZI 

LaFeOs 

Orthorhombic 

m 

LaCoOs 

Rhombohedral 

m 

LaNiOs 

Rhombohedral 

na 

LaCuOs 

Rhombohedral 

m 


TABLE I: Low temperature phases of the LaB 03 perovskites 
in the ground state. 


FIG. 2: (Color online) A molecular orbital diagram for an 
ABO3 perovskite. The bonding states (7r and a) have pre¬ 
dominant O 2 p character, the antibonding states (n* and cr*), 
have predominant B-cation d character. For simplicity, only 
the valence states are included, and the A-cation orbitals have 
been neglected. 


to a greater hybridization with O 2 p states, suggesting 
the 7r-bonding p — t2 g interaction is dominant in driving 
ferroelectricity; however, both the d e g states and the O 
2 p states are axial in nature, and therefore have a greater 
potential for overlap (and hence covalency), and therefore 
the cr-bonding p — e g interaction is also important. 


B. The LaBC>3 Series 


The key goal of this study is to determine the mi¬ 
croscopic origins of the ferroelectric and multiferroic be¬ 
haviour of BiCoOs; however, based on our previous work 
m , we propose that the Co 3+ ion drives the lattice insta¬ 
bility, rather than (or in conjunction with) the Bi 6 s hy¬ 
bridization. To isolate the contribution of the B-cation to 
the ferroelectric lattice instability in ABO 3 perovskites, 
we perform a trend study of the LaBOs series (ruling out 
contributions from Bi). A comparison between the con¬ 
tributions of Bi and La to the ferroelectric instability of 
perovskites has been previously reported [T| ; additionally, 
we will further compare the role of the A-cation in stabi¬ 
lizing the polar phase for LaBOs and BiBOs perovskites 
in section III C| For the remainder of this section, we fo¬ 
cus on the role of the B-cation in the lattice instability 
of LaBOs perovskites. 

As has been discussed above, perovskites will typically 
undergo a competing centrosymmetric distortion, rather 
than a polar one; this is true across the entire LaBOs 
series. In Table [l| we present the experimentally re¬ 
ported ground state phases for each LaBOs perovskite. 
As can be seen, the LaBOs perovskites typically prefer to 
undergo orthorhombic and rhombohedral distortions in¬ 


volving the tilting and rotation of B 06 octahedra, rather 
than exhibiting polar ferroelectric distortions. This fact, 
and the general lack of magnetic ferroelectric perovskites 
PQ, necessitates the study of artificial systems so as to 
perform trend studies that can allow insight into the fer¬ 
roelectric d° rule, and perovskite multiferroics. Histor¬ 
ically, one of the key ways in which multiferroicity has 
been studied, is to take a magnetic perovskite like the 
LaBOs systems presented in Table [T| and to calculate the 
ferroelectric lattice instability in the cubic phase, e.g., by 
analysis of phonon softening mm- For example, anal¬ 
ysis of the spin - phonon coupling in the cubic phase 
of LaMnOs and LaFeOs has been used to reveal multi¬ 
ferroic effects [32]. With this same approach, below, we 
perform a qualitative trend study to examine the ferro¬ 
electric lattice instability in the LaBOs perovskites as a 
function of the d orbital occupation and configuration, 
so as to provide insight into perovskite ferroelectricity 
and possible multiferroicity; for simplicity, to begin with, 
we neglect the long range inter-ionic magnetic ordering 
for the LaBOs systems as we try to elucidate qualitative 
trends. However, in the next section, we present quanti¬ 
tative calculations focusing on magnetoelectric coupling 
for the realistic perovskite system BiCoOs. 

In this section we study the LaBOs series, where B 
represents all of the possible B-cations (in the 3+ charge 
state [42]) across the 3d-series (i.e., B = Sc 3 + (d°), ..., 
Cu 3 + (d 8 )). The electron configuration for each B-cation 
under study is presented in Table [ h| For d A —d 7 B-cations, 
the octahedral crystal field splitting leads to two possible 
orbital configurations, known as the low spin (LS) and 
high spin (HS) states, which maximize and minimize spin 
pairing, respectively (in some cases an intermediate spin 
(IS) state is also possible). As representative examples 
of the electronic structure for the LaBOs systems, the 
density of states (DOS) is plotted for LS and HS LaFeOs 
in Fig. [3] as well as the expected orbital filling for the Fe 
3d band. As can be seen, the DOS confirms the orbital 
fillings of Table |iT| For LS LaFeOs, the d 5 configuration 
is such that only the t 2 g band is occupied, whereas the e g 
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band is empty; on the other hand, in the HS state each 
t 2 g and e g orbital is singly occupied, and there exists 
a strong exchange-splitting, as well as increased p — d 
hybridization. 

As previously discussed, in the simple cubic perovskite 
structure, the B-cation is 6-fold coordinated to oxygen, in 
a configuration that is symmetrical with respect to inver¬ 
sion. We study the tendency for the LaBOs perovskite 
to exhibit a ferroelectric distortion, by calculating the 
zone-center phonon modes for this high symmetry cubic 
phase - a positive frequency indicates stability against a 
spontaneous distortion, an imaginary frequency, or soft 
mode, indicates that the system is unstable [43]. Upon 
calculating the eigenvectors of the dynamical matrix, we 
identify the zone-center phonon mode that has been pre¬ 
viously associated with ferroelectricity - it is the soften¬ 
ing of this mode that indicates a ferroelectric instability 
[I]. An example of the eigenvector for this mode, calcu¬ 
lated for the HS state of LaFeCU is presented in Table [111} 
our calculation closely reproduces the eigenvectors pre¬ 
viously reported for LaBOs perovskites [I]. In Table [TT} 
the frequencies of the ferroelectric mode (ufe) are pre¬ 
sented for each LaBOs system. As a further test, the lat¬ 
tice vectors and internal coordinates are presented for the 
LaBOs systems under study. The relaxations are started 
from a nearly cubic structure, however a small ferroelec¬ 
tric distortion is present. The systems that favour a fer¬ 
roelectric distortion relax so as to increase the internal 
polarization and the c/a axis ratio. On the other hand, 
systems which do not favour a ferroelectric distortion re¬ 
lax back to a paraelectric (cubic) structure, without in¬ 
ternal distortions and with c/a = 1. We point out that 
for an orbitally degenerate configuration of the d-band, 
the system is also unstable against a Jahn-Teller driven 
distortion [22]; however, this is outside of the scope of 
the current work, and we present results for the cubic 
structure of LaBOs systems that cannot be stabilized in 
the PJT-driven ferroelectric phase. 

Looking at Table [nj starting with LaScOs, a strong fer¬ 
roelectric instability is present - the ferroelectric phonon 
mode is softened (c ofe = 155.7i cm -1 ), and the sys¬ 
tem is strongly tetragonal (c/a = 1.25) with ferroelec¬ 
tric distortions of the Sc06 octahedron. This result is 
somewhat expected as Sc 3+ has the typical d° configura¬ 
tion for a ferroelectric B-cation. Moving across the 3d- 
series for Ti 3+ , V 3+ , and Cr 3+ , as electrons are added to 
the B-cation, the ferroelectric instability begins to dis¬ 
appear, and as the ferroelectric mode is hardened the 
paraelectric structure is preferred. The observed trend 
with increasing occupation of the d band is entirely as 
expected, and in accordance with the d° rule. What is 
somewhat surprising is that for HS d 5 —d 7 and d 8 cations, 
the ferroelectric instability is recovered. For HS LaFeCU, 
LaCoCU, and LaNiCU, and for LaCuCU, the ferroelectric 
mode is softened, and the polar distortion becomes en¬ 
ergetically favoured over the paraelectric phase. The HS 
state of LaMnCU exhibits a slight ferroelectric distortion 
(c/a = 1.02), with a total energy only slightly lower than 



FIG. 3: (Color online) In (a), the projected partial density of 
states is plotted for low spin (LS) (top) and high spin (HS) 
(bottom) LaFeCU The dotted red line represents O 2 p states, 
the solid blue and purple lines represent respectively the t 2 g 
(d xy , d xz , d yz ) and e g (d z 2 , d x 2 _ y 2 ) orbitals of the Fe 3d band. 
In (b) and (c), the d orbital configurations are shown for LS 
and HS Fe 3+ . 


the paraelectric phase; as described in section [HI A[ this 
system can be thought of as type (ii), with Ko ~ K v . 

It is noted here that the apparent coupling of the elec¬ 
trical polarization to the spin state is enormous - for the 
d 5 — d 7 cations, switching between spin states leads to 
a phase change between a paraelectric and a ferroelec¬ 
tric structure, i.e., these cations exhibit the multiferroic 
crossover effect [8]. For example, in Fig. [4] the crystal 
structure of LS and HS LaFeCU is shown (as obtained 
in Table 0 - In the LS state, the system is cubic, in the 
HS state, the unit cell is tetragonal and the atomic posi¬ 
tions show strong polar distortions. To further illustrate 
this, also in Fig. [4} the ferroelectric nature of LS and HS 
LaFeCU is demonstrated by moving the Fe ion from the 
central (paraelectric) position to an off-center (ferroelec¬ 
tric) position. The LS state of LaFeCU exhibits a single 
well profile with a global energy minimum at the central 
position. On the other hand, the HS state exhibits a deep 
double-well profile with minima at the off-centred posi¬ 
tions, similar to what is displayed by other ferroelectric 
materials such as BaTiOs and PbTiOs [4]. At this point, 
it is somewhat counter-intuitive to think that the ferro¬ 
electric lattice instability is enhanced in the HS state of 
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B-cation 

Spin state 

Elec. Config. 

ojfe (cm 1 ) 

a (A) 

C (A) 

La 

B 

01 

02 

A Efe (meV) 

Sc 3 + 


d° 

155.7z 

3.893 

4.783 

0.000 

0.423 

0.834 

0.294 

-327 

Ti 3 + 


d 1 ', t2g( t) 1 

46.6z 

3.950 

3.959 

0.000 

0.456 

0.952 

0.437 

-9 

v 3 + 


d 2 ; t 2g { t) 2 

3.0 i 

3.863 

4.023 

0.000 

0.484 

0.983 

0.474 

-6 

Cr 3+ 


d 3 ; t 2g (t) 3 

53.9 

3.785 

3.785 

0.000 

0.500 

0.000 

0.500 

- 

Mn 3+ 

LS 

d 4 ; t 2g (-t) 3 t 2g (l) 1 

87.9 

3.809 

3.809 

0.000 

0.500 

0.000 

0.500 

- 


HS 

d 4 ; 

50.7 

3.880 

3.940 

0.000 

0.489 

0.988 

0.485 

-9 

Fe 3+ 

LS 

d 5 ; t 2g (V 3 t 2g (l) 2 

142.6 

3.786 

3.786 

0.000 

0.500 

0.000 

0.500 

- 


IS 

d 5 ; < 2 9 (t) 3 i 2 9 4 ) 1 e 9 (t ) 1 

80.9 

3.839 

3.839 

0.000 

0.500 

0.000 

0.500 

- 


HS 

d 5 ; i 29 (t) 3 e 9 (t) 2 

87.8 i 

3.800 

4.620 

0.000 

0.440 

0.842 

0.319 

-76 

Co 3+ 

LS 

d 6 ; t 2g tf) 3 t 2g (l) 3 

131.4 

3.775 

3.775 

0.000 

0.500 

0.000 

0.500 

- 


IS 

d 6 ; t 2g (t) 3 t 2g (4-) 2 e 9 (t ) 1 

97.7 

3.815 

3.815 

0.000 

0.500 

0.000 

0.500 

- 


HS 

d 6 ; < 2 9 (t) 3 i 2 9 4 ) 1 e 9 (t ) 2 

118.lz 

3.807 

4.577 

0.000 

0.443 

0.827 

0.309 

-79 

Ni 3+ 

LS 

d 7 ; <2 9 (t) 3 i2 9 4) 3 e 9 (t ) 1 

112.3 

3.815 

3.815 

0.000 

0.500 

0.000 

0.500 

- 


HS 

d 7 ; i 29 (t) 3 * 2 9 4 ) 2 e 9 (t ) 2 

47.4z 

3.758 

4.632 

0.000 

0.453 

0.841 

0.318 

-208 

Cu 3 + 


d 8 ; i 29 (t) 3 *2 9 4)%(t ) 2 

61.3z 

3.830 

4.010 

0.000 

0.468 

0.938 

0.414 

-15 


TABLE II: The electronic configuration (Elec. Config.) of each B-cation studied is shown, along with the frequency of the 
ferroelectric phonon mode (ujfe) for each LaB 03 perovskite in the cubic phase (at T). An imaginary frequency suggests a 
ferroelectric instability, a positive frequency suggests the system is stable against a ferroelectric distortion. The relaxed lattice 
vectors and internal coordinates are also shown for the LaBOs systems. Internal coordinates are given in lattice units along the 
c direction, which is the direction that ions are displaced in a ferroelectric distortion. For systems that prefer the ferroelectric 
phase, the energy lowering (with respect to the paraelectric phase) associated with the distortion is also presented as A Efe- 



La 

Fe 

01 

02 

HS LaFeOs 

-0.60 

0.25 

0.28 

0.50 


TABLE III: An example of the eigenvector for the phonon 
mode associated with the ferroelectric lattice instability. The 
present calculation is for the high spin (HS) state of LaFe 03 . 


these B-cations, since here the local magnetic moments 
are maximal; in this case, it appears that a large num¬ 
ber of unpaired d electrons actually drives ferroelectricity, 
rather than inhibits it. 

With it now established which magnetic cations ex¬ 
hibit a strong tendency for a ferroelectric distortion, the 
next step is to develop a model as to why, in terms of the 
microscopic and quantum mechanical effects that drive 
magnetism and ferroelectricity. The systematic reduc¬ 
tion of the ferroelectric instability with orbital filling for 
the d° — d 4 B-cations is in accordance with the d° rule. 
A possible explanation for this effect has been discussed 
by Khomskii [2] in terms of two contributing factors: (i) 
the presence of unpaired d electrons will inhibit the for¬ 
mation of a p — d bonding spin singlet state (ti — it) 
due to the exchange interaction, and (ii) the population 
of the antibonding orbitals (with predominant d charac¬ 
ter) should also weaken the p — d covalent bond. How¬ 
ever, these explanations clearly break down for the case 
of HS d 5 — d 7 and d 8 B-cations. For example, HS Fe 3+ 


has 5 unpaired d electrons, and therefore the influence 
of the “pair-breaking effect” of the exchange interaction 
should be very strong, and moreover there is a significant 
population of the antibonding states, yet LaFeOs shows 
a strong ferroelectric instability. 

We identify that the stabilization of the polar phase for 
HS d 5 — d 7 and d 8 B-cations can be understood within 
the PJT theory for ferroelectricity; here, note that we 
refer to “stabilization” with respect to the cubic phase, 
not the ground state structures presented in Table [TJ Ac¬ 
cording to Eq. [3j the lattice is stabilized against a ferro¬ 
electric distortion via the Kq term. However, this term 
is expected to be strongly affected by an increase in the 
volume of the unit cell [44]. Indeed, looking at the lattice 
constants presented in Table [TT] for the HS states of the 
d 5 — d 7 and d 8 LaBOs perovskites the volume is signifi¬ 
cantly increased in the HS state. For example, LaFeOs 
undergoes a 23 % increase in the volume in the HS state, 
when compared to the LS state. This dramatic increase 
in the volume is due to the population of the antibonding 
(<r*) e g states. Recently, Garcia et al. [44] have shown 
that the local forces that stabilize against a PJT driven 
ferroelectric distortion are strongly reduced by increasing 
the volume. As an example, to determine quantitatively 
the role of the volume increase in stabilizing the polar 
phase of the HS state in the LaBOs systems, the fre¬ 
quency of the ferroelectric phonon mode ujfe for the LS 
state of LaFeOs has been calculated, however the volume 
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FIG. 4: (Color online) Crystal structure of LaFe 03 (LFO) in 
the low spin (a) and high spin (b) states. In (c), the total 
energy of LaFe 03 is calculated as a function of the Fe ion 
position in the unit cell, presented as a fraction of the c-axis 
lattice vector. Sr atoms are represented by green spheres, Fe 
are brown, O are red. 


has been fixed to the equilibrium volume of HS LaFeOs. 
Indeed, at the HS volume, the mode is softened signifi¬ 
cantly with oj fe = 72. li. This does indeed suggest a sig¬ 
nificant decrease in Ko for the increased volume. Keep¬ 
ing the same volume, and switching to the HS state, the 
mode is softened slightly further, with oofe — 87. Si. The 
discrepancy in oopp could be attributed to an increase 
in the contribution from K v . To further investigate the 
effect of the volume increase, we also performed relax¬ 
ations for the LS state fixed at the volume of the HS 
state; indeed, now a polar phase is preferred with re¬ 
spect to the cubic structure. The c/a ratio is 1.17 and 
the system exhibits internal ferroelectric distortions, fur¬ 
ther supporting the idea that a reduction in the Kq term 
stabilizes the polar state. The c/a ratio and the internal 
distortions are slightly less than those of the HS state for 
the same volume, which could also point to a contribu¬ 
tion from K v \ however, it seems that the effect of volume 
on Kq is particularly strong. 


With regard to the K v term, as discussed in sec¬ 
tion III from Eq. [5j it is clear that the strength of 
the vibronic coupling between a p and d orbital depends 
on several factors: (i) whether or not the promotion of 
an electron from the occupied p state into the unoccu¬ 
pied (or partially occupied) d state conserves the spin 


of the system [8], (ii) the matrix element, (p\5H/5q\d), 
which describes whether or not the distortion increases 
the p — d overlap, and (iii) the energy gap between the 
two orbitals ( E p — Ed ). It is this dependence that leads to 
the so-called multiferroic crossover effect, whereby differ¬ 
ent orbital configurations lead to an enhanced contribu¬ 
tion from the vibronic term [8]. Looking at HS LaFeOs, 
the ground state To has an orbital configuration with 
S = 5/2; an example excitation T n involving the e g band, 
that conserves the spin of the system, is shown in Fig. [5) 
In this case, the excited state has the same spin as To, 
and so contributes to the summation in Eq. [5j It follows, 
that the number of T n terms that contribute to the sum¬ 
mation in Eq. [5] (and the strength of the contribution) is 
dependent on the orbital configuration, i.e., the multifer¬ 
roic crossover effect 0, and therefore the p — d vibronic 
interaction (K v ) is enhanced in one spin state with re¬ 
spect to another. 

To confirm that the cr-bonding p — e g interaction is 
enhanced in the HS state, in Fig. [6] we plot the charge 
density isosurface for the e g states of the Fe 3+ ion in the 
LS and HS states of LaFeOs. In the LS state, the charge 
density for the e g molecular orbitals are localized on the 
Fe 3+ ion, without significant hybridization with the sur¬ 
rounding O 2 p states. On the other hand, in the HS 
state, the e g band is hybridized and has significant O 2 p 
character; this is a clear indication of an increase in the 
PJT p — d covalent interaction. A similar analysis shows 
an increased in the p — t^g hybridization. We are able to 
approximate the total increase in covalent bonding based 
on the T-point projections of the bonding valence states 
onto the Fe and O atoms. In the LS state, the Fe 3d 
states contribute approximately 8% to the bonding states 
of the predominantly O 2 p valence band. When switch¬ 
ing to the HS state, indeed the bonding states are now 
comprised of 16% Fe 3d character, confirming increased 
p — d hybridization. Therefore, to summarise, we propose 
that the recovery of a strong ferroelectric instability for 
HS d 5 — d 7 and d 8 B-cations is due to a significant reduc¬ 
tion of the Kq term due to the volume increase, as well 
as increased p — d bonding (K v ). 

As an additional point, we notice that the HS d 5 — d 7 
and d 8 B-cations cations all lie to the right of the 3d- 
block. LaScOs has a p — 3d band gap of nearly 5 eV 
(Fig# however moving to the right of the d-block, the 
increased nuclear charge should bind the d-electrons more 
tightly, and reduce the p — d energy gap; this would re¬ 
duce the denominator in Eq.[5j and increase the p — d 
vibronic interaction. In Fig. [5j we plot the energy dif¬ 
ference between the band center of the O 2 p band and 
that of the B-cation 3d band for the LaBOs systems [45] . 
As can be seen, to the left of the 3d-block the average 
energy of the B-cation 3d band lies at a much higher en¬ 
ergy than the O 2 p band. However, moving to the right 
of the 3d-block the energy gap is systematically reduced, 
and for d 5 — d 7 and d 8 cations the O 2 p and B-cation 
3d bands lie at a similar energy such that the potential 
for PJT driven ferroelectricity is greatly increased. How- 
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FIG. 5: (Color online) (a) On the left the ground state orbital 
configuration (To) is shown for high spin Fe 3+ in LaFe 03 . On 
the right, the same orbital configuration is shown however an 
electron has been promoted from the O 2 p band into the e g 
band of Fe 3+ , the label indicates that this is the lowest 
energy excited state involving promotion of an electron into 
the e g band. As can be seen T 0 and have the same spin 
S = 5/2, and so T^ 9 contributes to the summation in Eq. [HJ 
(b) The energy difference between the O 2 p and B-cation 3 d 
band centres is plotted for the LaB 03 series. The position 
of the band center is approximated [45], and this plot should 
only be considered as a qualitative guide. 


LS LaFe0 3 




HS LaFeCF 


FIG. 6: (Color online) The charge density isosurface is shown 
for the e g band of low spin (LS) (left) and high spin (HS) 
(right) LaFeOs. The isosurface is set to 10% of the maximum 
value. La ions are prepresented by green spheres, Fe are brown 
and O are red. 


ever, we point out that, clearly, a small p — d energy gap 
alone is insufficient to drive a ferroelectric distortion, as 
the condition K v > K 0 must be satisfied. 

The above discussion for the LaBOs perovskites should 
be quite generalizable for these B-cations if the A-cation 
is replaced. To confirm this, we studied the ferroelectric 
properties of LS and HS CeFe 03 . This allows us to con¬ 
firm the key findings from above, namely: (i) a strong fer¬ 


roelectric instability for a half-filled e g band (which is ab¬ 
sent for other orbital configurations), and (ii) the strong 
coupling of the ferroelectric polarization to the spin state. 
Indeed, we find that for HS CeFeOs the ferroelectric 
phonon mode is softened with ujfe = 97.5i cm -1 , and 
the relaxed structure is tetragonal with c/a = 1.24 and 
strong polar distortions; for LS CeFeOs the ferroelectric 
mode is hardened to ujfe = 112.5 cm -1 and the system 
relaxes back the paraelectric phase. Therefore, the be¬ 
haviour of Fe 3+ in LaFeOs is well reproduced in CeFeOs. 
We emphasise that, again, CeFeOs does not crystallize 
in the tetragonal ferroelectric structure; rather, this sys¬ 
tem prefers an orthorhombic distortion [46]. However, 
the previous result illustrates that the ferroelectric and 
multiferroic properties presented for LaFeOs arise due to 
the d 5 Fe 3+ B-cation, and not some other effect as a result 
of the A-cation. 

At this point it has been demonstrated that LaFeOs, 
LaCoOs and LaNiOs can potentially exhibit the multifer¬ 
roic crossover effect, and this strong coupling of the elec¬ 
trical and magnetic polarization is highly desirable from a 
multiferroic device perspective. However, unfortunately, 
these materials prefer to undergo centrosymmetric dis¬ 
tortions, rather than a polar one. Specifically, at room 
temperature LaCoOs and LaNiOs exhibit rhombohedral 
distortions [47, 48], whereas LaFeOs is orthorhombic [38] . 
The importance of competitive distortions in the d° rule 
has already been stressed previously [T]. Such distortions 
can originate from the JT effect, or, in the absence of a JT 
distortion, the PJT effect can also give rise to centrosym¬ 
metric distortions [49 j . One possible route to suppressing 
these competing distortions is via strain. Tensile strain 
is known to inhibit orthorhombic and rhombohedral dis¬ 
tortions, and favours a polar tetragonal phase [50]. It is 
therefore possible that through substrate engineering of 
thin films, the LaBOs systems studied above could be 
stabilized in the polar tetragonal phase via strain. In¬ 
deed, it was shown that strain engineering could be used 
to stabilize a weak ferroelectric state in magnetic AMnOs 
perovskites GEO Eg. a second approach to suppress¬ 
ing the competing centrosymmetric distortions is via the 
choice of the A-cation. It is a well known empirical fact 
that the tilting and rotation of B06 octahedra depends 
strongly on the relative sizes of the A and B cations 1531 • 
A larger A-cation compared to a smaller B-cation will 
inhibit these distortions. For perovskites of the A 3+ B 3+ 
charge ordering, La 3+ is already very large, and replacing 
the A-cation with most of the alternative 3+ cations (e.g., 
other the rare earth ions) would increase the tendency for 
competing centrosymmetric distortions. One possibility 
is to focus on A 2+ B 4+ Os systems, as 2+ cations are ex¬ 
pected to be larger - indeed, the well known ferroelec¬ 
tric perovskites PbTiOs and BaTiOs have the A 2+ B 4+ 
charge ordering. Surprisingly, there has been far less 
work on magnetism, or multiferroicity, in the A 2+ B 4+ Os 
perovskites. Focusing on d -fillings which can lead to the 
multiferroic crossover effect, i.e., d 5 — d 7 , and taking into 
account that Cu will not form in the 4+ charge state, the 
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FIG. 7: (Color online) The total energy (per formula unit) 
of cubic high spin (HS) BiCo 03 and cubic HS LaCo 03 is 
plotted as a function of the A-cation displacement from the 
high symmetry site. 

possible B-site cations are Co 4+ and Ni 4+ . For large A- 
cations in the 2+ charge state Sr 2+ , Ba 2+ and Pb 2+ are 
plausible candidates. A survey of the literature suggests 
that in most cases the different combinations of these 
A- and B-cations will result in a ground state structure 
in the hexagonal phase [54H56] , which commonly occurs 
when an A-cation is much larger than the B-cation. For 
the case of SrCoOs, it was indeed reported that the cen- 
trosymmetric distortions were suppressed, and the cu¬ 
bic structure is the ground state phase 1523- However, 
SrCoOs has an IS ground state with a orbital fill¬ 
ing which, as we have shown, does not lead to a strong 
ferroelectric instability. 


C. Multiferroic Crossover: BiCoOs 

In the previous section of the manuscript, it was re¬ 
vealed that certain d n cations can indeed provide the 
driving force for a ferroelectric distortion; however, com¬ 
petition with noncentrosymmetric distortions can ac¬ 
count for the empirical d° rule [l, 5]. The perovskite 
BiCoOs appears to be an exception to this rule, as the 
ground state crystal structure is that of a tetragonal fer¬ 
roelectric phase, and yet the d 6 Co 3+ ion is in the HS 
state with (7-type antiferromagnetic ordering US] (for 
more detailed information on this system see Ref. [T6]). 
It was suggested by Jia et al. that hybridization of the 
high lying Bi 6s state with the O 2 p valence band stabi¬ 
lizes the ferroelectric distortion m • Indeed, surely the 
Bi ion will play a role. To demonstrate this, in Fig.[7| we 
plot the total energy of cubic HS BiCoOs as a function of 
the displacement of the Bi ion from the high symmetry 
position. For comparison, we also plot the same result 
for cubic HS LaCo 03 . Clearly, while the total energy of 
LaCoOs is increased by the displacement, the total en- 



FIG. 8: (Color online) The partial density of states (DOS) 
for low spin (LS) (top) and high spin (HS) (bottom) BiCoOs - 
The LS state of BiCoOs is non-magnetic, whereas the HS 
state has C-type antiferromagnetic ordering. The DOS has 
been shifted so that the Fermi level is at zero. 


Spin state 

a 

c 

Bi 

Co 

01 

02 

LS 

3.731 

4.056 

0.000 

0.426 

0.881 

0.360 

HS 

3.727 

4.825 

0.000 

0.424 

0.892 

0.269 

TABLE IV: 

Lattice 

vectors 

and 

internal 

coordinates of 


BiCoOs in the low spin (LS) and high spin (HS) states. 


ergy of BiCoOs is reduced, consistent with Bi 6s - O 2 p 
hybridization m • However, based on the results of Ta¬ 
ble [TTJ it seems likely that the HS state of a d 6 ion provides 
an additional driving force for the ferroelectric distortion 
via PJT vibronic coupling. Indeed, in a previous report, 
we have demonstrated the multiferroic crossover effect 
for Co 3+ dopants at the Ti site in PbTiOs [T7] . In the 
report by Jia et al. m , the authors fixed the volume of 
their BiCoOs cell to that of the experimental structure, 
and restricted the Co 3+ ions to be in the LS state; in this 
case the system favours a tetragonal phase, and based on 
this result, the authors concluded that the distortion is 
due purely to Bi - O hybridization, ruling out a contribu¬ 
tion from Co 3+ . However, clearly, by fixing the volume 
of the system to that of the experimental structure, the 
effect of the reduced Kq term caused by the LS - HS spin 
transition cannot be observed; additionally, the increase 
in the vibronic term for HS Co 3+ was not explored. 

To determine the role of the HS Co 3+ ion in driving the 
ferroelectric state of BiCoOs, we have investigated the 
atomic and electronic structure of BiCoOs using accu¬ 
rate hybrid functional calculations. Our results confirm 
that BiCoOs has a HS configuration as the ground state 
with (7-type antiferromagnetic ordering, and a strongly 
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FIG. 9: (Color online) The calculated crystal structure for 
BiCo 03 in the low spin (a) and high spin (b) states. In (c), the 
ground state C-type antiferromagnetic ordering of the high 
spin Co ions is shown. Bi atoms are represented by purple 
spheres, Co are blue and O are red. 


tetragonal {c/a = 1.29) crystal structure with large in¬ 
ternal ferroelectric distortions. Our calculated crystal 
structure for the LS and HS states of BiCoOs is shown 
in Fig. [9j the C-type antiferromagnetic ordering is also 
shown. A \[2a x a/26 x 2c, 20 -atom BiCoOs cell is used to 
account for the long range magnetic ordering. The lat¬ 
tice vectors and internal coordinates calculated for the 
LS and HS sates of BiCoOs are presented in Fig. [IV] Im¬ 
portantly, we find that in the LS state the axis ratio is 
reduced to c/a = 1.08, and the internal ferroelectric dis¬ 
tortions are reduced; this result confirms the multiferroic 
crossover effect for BiCoOs, and moreover, this result 
reveals the crucial role of the Co 3+ ion in driving the 
ferroelectric state, which has been previously attributed 
to Bi 6s - O 2 p hybridization. It is noted, however, that 
unlike the LaBOs systems, BiCoOs exhibits a remnant 
polarization in the LS state, supporting the idea that 
the Bi 6s does indeed provide a stabilizing effect for the 
ferroelectric distortion. 

As the ferroelectric polarization of BiCoOs is affected 
by switching the spin state, potentially, manipulation of 
the polarization - e.g., by strain, or by an externally ap¬ 
plied electric field - could induce spin crossover. This 
idea is supported by recent first-principles calculations 
based on the generalized gradient approximation (GGA) 
m, where the paraelectric phase of BiCoOs was re¬ 
ported to prefer the LS state. However, the LS state 
was found to be metallic which would prevent switch¬ 
ing via an applied electric field. It should be pointed out 
though, that these predictions for the electronic and mag¬ 
netic structure should not be considered reliable as the 
GGA cannot accurately predict the band gaps of semi¬ 
conductors, and because the GGA will overly favour the 
LS state due to the large self-interaction error [21]. 

Using advanced hybrid functional calculations, in 
Fig- m we are able to confirm that BiCoOs is semicon¬ 
ducting in both spin states, with band gaps of 2.3 eV and 
1.2 eV for HS and LS BiCoOs, respectively. This result 
suggests that an applied electric field could be used to 
manipulate the internal polarization. The hybrid func¬ 
tional calculations are also able to give an accurate treat¬ 



FIG. 10: (Color online) The total energy (per formula unit) 
of low spin (LS) and high spin (HS) BiCoOs is plotted as 
a function of the generalized coordinate. To model the de¬ 
pendence of the spin state energetics on the polarization, the 
internal coordinates are moved continuously from those of the 
relaxed high spin structure (generalized coordinate = 1) to a 
paraelectric structure (generalized coordinate = 0). 


ment of the HS - LS splitting of spin crossover systems 
EHE2], and therefore we are able to more accurately pre¬ 
dict the way in which changes in the electric polarization 
will affect the spin state energetics. An applied field can 
reduce the internal polarization of a ferroelectric mate¬ 
rial [58] , and force the internal coordinates to be closer to 
that of the paraelectric phase. In Fig. [To] the dependence 
of the spin state energetics on the ferroelectric polariza¬ 
tion is presented for BiCoOs. The internal coordinates 
are fixed by interpolating between that of the relaxed HS 
geometry and the paraelectric phase, however at each 
point we allow relaxation of the lattice vectors. As can 
be seen, at the polar HS geometry, the total energy of the 
HS state is far lower than that of the LS state. As the 
internal polarization is reduced, the HS - LS splitting is 
reduced and at about the halfway point of the interpo¬ 
lation, the LS state becomes the ground state (note that 
we also performed test calculations for other magnetic 
configurations of the HS state, and also for the IS state 
to confirm that the LS state becomes the ground state). 
Figure [lO] provides a good first indication for electric field 
control of magnetism in BiCoOs, i.e., changes in the in¬ 
ternal polarization strongly affect the relative energies of 
the HS and LS states; we propose experimental studies 
to confirm this magnetoelectric effect. 


IV. SUMMARY 

In this report, using the LaBOs series as a model, we 
have performed a qualitative trend study so as to inves¬ 
tigate the interaction between ferroelectricity and mag- 
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netism in ABO 3 perovskites. For the B-cations stud¬ 
ied, moving to the right of the 3d-series, it was found 
that initially, increasing the occupation of the B-cation d 
orbital decreases the tendency for a ferroelectric distor¬ 
tion, which is in agreement with the ferroelectric d° rule. 
However, a surprising result was found for HS d 5 — d 7 
and d 8 cations, in that a strong ferroelectric instability 
was recovered. This result was explained in terms of the 
pseudo Jahn-Teller theory for ferroelectricity, and it was 
demonstrated that, contrary to the current understand¬ 
ing of the ferroelectric d° rule, in some cases unpaired 


d electrons actually drive ferroelectricity, rather than in¬ 
hibit it. For the case of BiCoOs, the crucial role of the 
Co 3+ ion in driving the ferroelectric lattice instability has 
been revealed. Finally, it was demonstrated that d 5 — d 7 
B-cations will exhibit the multiferroic crossover effect, 
whereby switching between spin states can strongly af¬ 
fect the ferroelectric polarization; moreover, the manip¬ 
ulation of the polarization will also affect the relative 
energies of the different spin states, which suggests the 
possibility of electric field control of magnetism. 
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